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Abstract. This paper deals with the two-dimensional incompressible, laminar, steady-state boundary layer 
equations. First, we determine a family of velocity distributions outside the boundary layer such that these prob- 
lems may have similarity solutions. Then, We examine in detail new exact solutions, called Pseudo-similarity, 
where the external velocity varies inversely- linear with the distance in the x-direction along the surface {Ue{x) = 
UooX~^). The analysis shows that solutions exist only for a lateral suction. Here it is assumed that the flow 
is induced by a continuous permeable surface with the stretching velocity U^x^^. For specified conditions, we 
establish the existence of an infinite number of solutions, including monotonic solutions and solutions which 
oscillate an infinite number of times and tend to a certain limit. The properties of solutions depend on the suc- 
tion parameter. Furthermore, making use of the fourth-order Rungc-Kutta scheme together with the shooting 
method, numerical solutions are obtained. 
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1. Introduction 

In this paper we are concerned with the classical two-dimensional laminar incompressible boundary layer 
flow past a wedge or a flat plate [221 ■ For the first approximation, the model is described by the Prandlt 
equations or the boundary layer equations 

(1.1) udxU -\- vdyU = UedxUe + ydyyU, dxU + dyV~0, 

where (x, y) denote the usual orthogonal Cartesian coordinates parallel and normal to the boundary y — (the 
wall),u(a;, y), v{x, y) arc the x and y velocity components, respectively, and > is the kinematic viscosity. The 
function Ue = Ue{x) is a given external velocity flow (the main-stream velocity) which is assumed throughout 
the paper to be nonnegative and is such that u{x,y) tends to Ue{x) as y — > 00. Equations can be written 
in the form 

(1.2) dy^Pdly^P - dx^d^yyij = U^Ue + 1^ d^y y y^ , 

where ip is the well-known stream function defined by u = dy^p, v = —dx'4^- 

This equation with appropriate external velocity flow has been the main focus of studies of particular exact 
solutions. It is well known that to derive properties of solutions to a nonlinear partial differential equations 
we use a family of special solutions. They play an important role for describing the intermediate asymptotic 
behavior of classes of solutions of original problems with arbitrary initial data (see for example @],[2I|). ^ 
crucial step, in the analysis, is to get favorable conditions such that particular solutions exist. For the boundary 
layer problems, research on this subject has a long history, which dates to the pioneering works by Blasius [Hj 
and Falkncr and Skan Their investigations lead to solutions to (1.2) in the form 

(1-3) ^{x,y) = x^fiyx-^). 

* Corresponding author. E-mail; zakia.hammouch@u-picardic.fr 



LAMFA, No.9/Mars 2004, p. 1 



Therefore, if is a such solution, it is easily verified that, for 7^ 0, 
for all A > 0, and we immediately see that 

This means that a solution ip{x^y) for y fixed is similar to the solution Tp(xo,y) at a certain xq- This solution 
is called invariant or similarity solution and the function / is called the shape function or the dimensionless 
stream function. 

The main goal of identifying similarity solutions is to reduce the original problem to an ordinary differential 
equation which is easier to analyze. One says that the function (1.3) scales the partial differential equation (1.2) 
if the function / satisfies an ordinary differential equation called similarity boundary layer equation. We refer 
the reader to [21jIZ1,|H1j [22] 1 HH] a-nd the references therein. A relation of the type a = h{P) called scaling 
relation. 

In |Hj the Blasius equation 

(1.4) /"' + i//" = on (0,cx3) 

is derived from (1.1) and (1.3) where Ue is assumed to be a constant function. In the above equation the primes 
indicate differentiations with respect to the similarity variable t = yx~^^^. Here we have a = /? = i. The Blasius 
equation is also a particular case of the Falkner-Skan equation |13| 

(1.5) + on (0,cx3), 

where the external velocity Ue{x) — UooX™', {Uoo > 0). The exponents a, (5 are given by a = ■^^^^i P — ^—r^- 
Guided by the results of [H], 23 attention will be given, in Section 2, to identify a class of external velocities 
for equation (1.2) to possess solutions under the form (1.3) where a + (3 = 1. Additionally, we investigate, in 
Section 3, the similarity solutions to (1.2), where both of the external velocity and the stretching velocity of the 
permeable surface are assumed to vary as x~^. 

2. Similarity solutions 

As it is said in the introduction, this work deals with the similarity steady boundary-layer flow induced by 
an incompressible viscous fluid past a semi-infinite flat plate. The phenomenon is governed by the system 

(2.1) udxU + vdyU = I'dyyU + UedxUe, dxU + dyV^O. 
accompanied by the boundary condition 

(2.2) u{x,y)^Ue{x) as y ^ 00. 
Therefore the stream function satisfies 

(2.3) dy^dlyi, - dxi^dlyi> = UAUe + Vdlyy^b, 

and 

(2.4) lim dyili{x, y) = dyipix, 00) = Ue{x). 

y— »oo 

The main problems, arising in the study of similarity solutions, are related to the existence of the exponents 
a and (3 and to the rigorous study of the differential equation satisfied by the profile /, which is, in general, 
nonlinear. For equation (2.3), the classical approach for identifying a and P is the scaling and transformation 
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group The essential idea is to seek a and b such that if ip satisfies (2.3) the new function ip\{x,y) = 
X°-ip{X''x, Ay) is also a solution. This will be certainly possible if the external velocity field Ue is subject to the 
transformation. The parameters a and b may depend on a and p. Using (2.3), it is easily verified that tpx is a 
solution, for any A, if the following 

(2.5) a + b^l, X^+''Ue{X^x)Ue, (X^x) = Ue{x)Ue^ [x), 

hold for any x (see the proof of Theorem 2.1 below). In Falkner and Skan considered the case where 

(2.6) Ue{x) = U^x"-', 

where Uoo{> 0) and m are constants. The case m = was treated earlier by Blasius [Hj. Therefore, we deduce 
from (2.5), for to ^ 1, 

m + 1 ^ 2 

a= , 0= . 

TO — 1 m — 1 

Now assume that if; satisfies the invariance property ipx = ip. Thus if we set X^'x = 1, wc find that tp can be 
written in the form 

^p{x,y) = x°'f{yx~'^), 

where 

This leads to (1.3) with a ~ ^2i±i and /? = ^'^^^y^- This, may be, helps us to understand the similarity stream 
functions obtained by Blasius and Falkner and Skan. These authors considered solutions to (2. 3), (2. 4) of the 
form 

(2.7) i; = {iyUex)^f{t), t = {Ujyx)-^y. 

In fact, we shall sec in the following results, that condition (|2.()|) is necessary for problem (2. 3), (2. 4) to admit 
similarity solutions under the form (1.3). The main result of this section says that (1.3) scales (2.3) if the 
external velocity satisfies 

t/e(x) =C/oo(.T+2;o)", 

for some xq, where m = a — (3. We note, in passing that if ip is a solution the function 'il'{x, y) ~ ip{x + xq, y) 
satisfies (2.3) with Ue{x + xo) instead of Ue{x). Then it allows xq to be zero. The necessary condition on Ue 
can now be stated as follows. 

Theorem 2.1. Assume that equation (2.3) has a similarity solution in the form (1.3) where a + f3 = 1. Then, 
there exist two nonnegative constants, ci , C2 sucii that 

(2.8) C/f(x) =Cix2"+C2, 
for all X > 0, where m, = a — (}. 

Proof. Let %p he & stream function to (2.3) defined by (1.3) where a + /3 = 1. Assume first that /? 7^ 0. We 
choose a = — ^, b ~ ^, and define ip\{x, y) = X°'ip {X^x, Ay) . Hence a + b = l,t/j = ip\ and 

L{^^){x,y) = X-+^Lmx''x,Xy) 
for any A > 0, where L is the operator defined by 

According to equation (2.3) we deduce for any A > and any x > 0, 

h{x) = X''+^h{X^x), 
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where h(x) = Ue{x)dxU{x). In particular, for fixed > 



Setting x ~ X'^xq wc infer 
Solving the equation 

(2.9) = x-^xp^h{xo) 

yields us (2.8) for /3 7^ 0, since + 1 = 2(a - /?). 

For /3 = 0, hence a = ?n = 1, the new function 

Mx,y) = X-ij{X--x,y), 

for any fixed a 7^ 0, is equivalent to ijj and satisfies 

L{ijx){x,y) = X-LmX--x,y) 

for any A > 0. Arguing as in the case /3 7^ one arrives at (2.8) with jti = 1. ■ 
A similarity assumption to (2.9) was proposed by Spalding and studied extensively by Evans |12) . The 
authors assumed that Ue satisfies the equation 

(2.10) ^ = CC/2(''-i)/^ 

ax 

where C is a constant. The similarity equation is then 

(2.11) /"' + //" + 7(1 -/'') = 0, 

which is equivalent to (1.5) with 7 = Now let us discuss the consequence of the boundary condition at 

infinity. We note, in passing, that from H2.8|l and 12.10|) we easily deduce that Ue{x) = UaoX™'. 

Theorem 2.2. Equation (2.3) has a similarity solution, ip, in the form (1.3) such that a+f3 — 1 and dyip{-, 00) = 
Ue, if and only if there exists a real m such that Then 

(2.12) Ueix) = U^x"\ 

for all a; > 0, Uoo is a constant. Moreover the real m satisfies m ~ a — [3. 

Proof. Let Ue{x) = Uoox"^. the existence of similarity solutions in the form (1.3) where a = ^^^^ and 
/3 = is given in |S] and ^31 (see also ^21,^11) for some values of m. Conversely, assume that (2.3), (2.4) has 
a similarity solution in the form (1.3) with a + f3 = 1. Then Ue{x)^ = cix^™ + 02 and limy^oo a:^'"(/'(?/x^^)^ = 
cix^™ +C2. Thereafter, the function /' has a finite limit at infinity, which is unique and is given by ci +C2X~'^"^. 
This is acceptable only for C2 = 0. ■ 

Remark 2.1. It may be noted that in the above theorem it is not required to find the range of the real m such 
that problem (2. 3), (2. 4) has a similarity solution in the form (1.3). In fact, the similarity equation may have 
no solution for some real m (see Section 3). Theorem 2.2 indicates, in particular, that for a prescribed external 
velocity Ue{x) = Uoox"^, the reals a and /3 such that (1.3) scales (2.3) are given by 

m + 1 TO — 1 
" = — 7; — 'P = — • 
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Remark 2.2. For a general external velocity, our approach can be used to obtain particular solutions. In ^5] 
it is obtained particular solutions, having the form 

tp ^ H{x, y)f{-q{x, y)) + Vo- 

So, if 

Ue{x)=Ax'''+U{x) 
we conjecture that there exist solutions in the form 

il^ix^y) = x°'f{yx~'^) +ijo{x,y), 

where a = ^^^^ and (3 = — ^^^j^- The function tpo may be similarity and connected to x" f{yx~^) and U{x). In 
the next result we shall identify the external velocity such that equation (2.3) has a solution in the form 

(2.13) V(a;,y) =a;"/(ya;"'^)+cx"y, « + 1, const. 

Note that the function ipoi^, y) = cx"'y can be written in the form "00(3^, y) = cx"^ {yx~^^), for any ai, /?i such 
that ai — /?! = n; that is i/'o is similarity. So, for a — P = n (2.13) is identically to (1.3). Therefore we assume 
that a ~ [3 ^ n. 

Theorem 2.3. Let a, (3 axe reals such that a — (3 ^ n. The reals a, (3 in (2.13) scale (2.3) if and only if 
a=|,/3=-|, n = — i and the following 

Ue{x) = Ax^ + cx~^, 

holds for some constant A. 

Proof. For Ueix) ~ ax^^^ + cx~^^^, it is shown in |2S] that the function ?/'(.t, y) = x"^/^ f{yx~^/^) + cyx~^/^ 
satisfies (2.3) 

Now, assume that (2.3) has a solution in the form (2.13). Substituting (2.13) into (2.3) yields 



(2.14) 



vf" + aff" -{a-l3)f''^ + c(7i + a-/3)a;"+"-'3-i/' 

+ncln + /?)ta;"+"-^-i/" - cnx'^'^-^ + U^'^^^ = 0, 



where t = yx ^ . Since a — (3 ^ n, the function / satisfies an ordinary differential equation if and only if 
(2.15) n + a- [3 = n + (3 = 0, 

and there exists a real 7 such that 

+ «//" - (a - /?)/'' = 7, 
From (2.13) and (2.14)2 we deduce immediately that a = = i, n= — -i and Ueix) = Ax^ + cx^^ . 



(2.16) 



Remark 2.3. Let us now derive the well known Blasius and Falkner-Skan. Of course the external velocity is 
given by Ue{x) = Uoox"^- We recall that for the Blasius model we have m = and the case m ^ was considered 
by Falkner-Skan. Instead of taking a = ^^^^ and /3 = ^^^^5-^ we shall insert (1.3) into (2.3) and choose a and 
f3 such that / satisfies an ordinary differential equation. Obviously we shall obtain that a and f3 must to be , 
respectively and —^^^y^. 
Problem (2. 3), (2. 4) is written as 

(2.17) dy^pd^^P - d^i^dly^ = vdlyyi, + inUlx^^-\ 
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with the boundary conditions 

(2.18) d^ipix,0) ^ dyip{x,0) ^0, dytP{x,oo) = U^x"" . 

If we substitute (1.3) into the first equation of this problem we obtain 

(a - /3)x2(«-/3)-i(j/)2 _ ax^{c^-fJ)-iff _ a;"-3/3^/// _ ^^2^^2 

and this is an ordinary differential equation if and only if 

2(a - /3) - 1 = a - 3/3 = 2m - 1, 



(the scaling relation) i.e. 



m + 1 , ^ m — 1 

and p = . 

2 ^ 2 



After a scaling ( ^/vUoof ^- y ^ j instead of /(.)) the corresponding ordinary differential equation is 

(2.19) f" + !!l±lff" + mil -f'') = 0. 

'u, 



The prime denotes the derivative with respect to i = \/^^yx 2 . The boundary conditions read 
(2.20) /(0) = /'(0) = 0,/'((^) = l. 

Remark 2.4. We observe, in passing, that the boundary conditions on the plate are not required. This means 
that our analysis works even if we have a continuous stretching surface. In this case, if the stretching velocity 
is given by Us{x) = U^x^ , we deduce from (1.3), n = a — (3 which is the same exponent as that of the external 
velocity. If u = at y = the function / satisfies |17 | .|18 | 

/"' + ^//"-rn(/''-l) = 0, 

/(O) = 0, /'(O) = C, /'(O) = 1, 
where C is equal to the ratio of the free stream velocity to the boundary velocity. 

In the case where the external velocity is zero, the stretching and the suction/injection velocity have the form 

u{x,0) = Uwx"\ v{x,0) ^ ~VujX^^ , 

the similarity equation is given by 



(2.21) 



/"' + 2i±i././" - m/'^ = 0, 
/(O) = a, /'(O) = 1, fi^) = 0, 



where a is a real (the suction/injection parameter). This problem arises also in the study the free convection, 
along a vertical flat plate embedded in a porous medium. We refer the reader to the papers [2^1 1 ISl i El i CHI and 
the references therein for a complete physical derivation and analysis of this problem. 

Remark 2.5. We finally, mention a result of 130] which provides that for the external velocity satisfying (2.12) 
where Uoo > and m > 0, any solution u to (1.1) with = Vw = such that Uy is continuous for < a; < 
00, < ?/ < 00, satisfies the fundamental asymptotic estimate 



1 + m log X 



as a; ^ 00, uniformly in y, where / is a solution to (2.19);(2.20). 
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3. The pseudo-similarity solutions 



In the present section we focuss our attention to the case m = — 1 . The case m 7^ — 1 has been abundantly 
studied. In particular it was considered for the Falkner Skan equation (2.11) where 7 = 2m/(m+ 1). The cases 
7 = and 7=1/2 are referred as Blasius and Homann differential equations, respectively. In Coppel 
classified all solutions of H2.11(l . where < 7 < 2. Craven and Pclctier ^0! proved that equation (2.11), where 
i < 7 < 1, has at most one solution / satisfying the boundary condition 

(3.1) /(O) = /'(O) = 0, /'(c^) = 1. 

In Hasting and Siegel showed that for I7I sufficiently small there exists a unique solution to (2. 11), (3.1) 
satisfying |/'| < 1 and /"(O) < 0. The case I7I > 1 can be found in the works ^| ^] by Hasting and Troy. 
In J2| deals with 7 > 1. In this it is shown that the equation has a periodic solution and for any integer N 
problem (2. 11), (3.1) has a solution with at least N relative minima. It is also mentioned that as 7 increases, 
the structure of periodic solutions and solutions to (2. 11), (3.1) gets progressively more comphcated. 
In the present section we restrict our attention to the case rn = —1. Therefore we shall consider the external 
velocity is given by Ue{x) = ^^(m = — 1), and get new solutions to the problem 

subject to the boundary conditions 

(3.3) dy^{x,0) = UwX~^ , dxi^{x,0) ^ —VioX~^ , dyij} {x, co) = U 00 ■x~^ , 

where Vw is a real and Uw and Uoo are nonnegative and satisfy Uw < Uoo- The subject of the present section is 
motivated by the work of Magyari, Pop and Keller j^T], where the case Uoo = was considered and discussed 
in detail. The authors showed that if m = —1 problem (2.21) has no solution. In order to overcome this 
difficulty the authors showed that the term Vu,log{x) must be added to expression (1.3). It has been also 
confirmed, by numerical calculations, that new solutions (pseudo-similarity solutions) exist provided that the 
suction parameter is large. 

For (3. 2), (3. 3), where Uoo 7^ and according to Section 2, the function ip can be written as 



ipix,y) = y/uUZfi^J^^yx ^). 



Since, in general Vu, = ^-^^^^ \/ vU 00 f {^) , we deduce = for m = — 1 and the function 9 = f satisfies 

-1 = 0, 



2 

3" I a'2 



0(O) = C,e(oo) = l, 



(3.4) 

where C = jj^ is in the interval [0, 1). 

The stability of equilibrium point (1,0) of (3.4) cannot be determined from the linearization. To analyze the 
behavior of the nonlinear equation (3. 4)1, wc observe that 

E'it) = 0, 

where E is the Lyapunov function defined by 
Then, for some constant C, the following 



±V2 ( C 



1 \ 1/2 



holds. The analysis of the algebraic equation of the phase path in the phase plane {0, 9') reveals that the 
equilibrium point (1, 0) is a center. Hence Problem (3.4) has no solution for any ^ > —1 except the trivial one 
61 = 1, (see Fig.3.1). 
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Fig. 3.1 Classification of solutions oi 9" + 6^ - 1 = according to 61(0) and e'(0). 



Remark 3.1. If we impose the condition 9{oo) = —1 instead of 9{oo) = 1-which is also of physical interest- it 
is easy to see that for any C ^ 2 there exists a unique solution-up to translation-. This solution satisfies 

and we find that 

6'(t) = 2- 3tanh^ ±i/\/2 + arctanh|(2 - C)/3)^/^} . 
To obtain exact solutions to (3. 2), (3. 3), we introduce the following Anzats for the stream-function 
(3.5) ij{x, y) = ax"F{x, byx'^). 

Guided by the analysis of Section 2 we take a = ^^^^ ^^^^ ■, o- = \fvU^ and h — . Here the real m is 

assumed to be any real. Hence boundary condition (3.3) reads 

dy^(x,^) = U.^„x^ , d^i}{x,^) = -V^,x^ , dyip{x, oo) = Uoox"\ 
Using this and the equation of the stream function we find that F satisfies 

F'" + ^FF" - m{F'^ - 1) + X {F"d^F - F'd.^F') = 0, 



(3.6) 



F'{x,0)=C, F'{x,^) = l, 
^^'^-F{x,0) + xdF{x,0) ^ 



where the primes denote partial differential with respect to t = \ ^^yx^~ . By writing 



F{x,t)= f{t)+H{x), 



we find 



(3.7) 



/"' + - m - l) + /" {xH' + i±™i/) = 0, 



'-±^m+^A^H{x)^xH'{x) 



=, /'(0) = ce [0,1), /'(oo) = i. 
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Hence, there exists a real t such that 

' /"' + H^//" -m(r-l]+ rf" = 0, t>0, 

(3.8) 



xH' + ^P^H = r, X > 0, 



2 

i±^/(0) + r=-^, /'(0) = C, /'(oo) = l. 
First, let us note that (2.19) and (3.8)i are equivalent for m ^ —1. In this case the general solution of (8.8)2 is 



1 + m ' 



where C is a constant. Thus the stream function is given by ip{x,y) = aC + ax ^ ( /(i) + ) i ^-i^d the new 

function g = f + satisfies equation (2.19). 
Thereafter, we will assume that m = — 1 and this leads to 



(3.9) 



r+Tf' + f' -1 = 0, 

/'(0) = C, /'(oo) = l, 



and 

H(x) = rlogx + C, C = const. 
Then the required exact solution has the form 

(3.10) ■4'{x,y) = af{hyx^^) + ar log(a:;). 

This formula is similar to the one studied in |24| in the context of rough surface growth and can be regarded as 
solution with dynamic scaling In passing, we note that from (3.8) one sees r = {vUoo) ■ This means 
that the constant r plays the role of suction/injection parameter. 

To study (3.9) it is more convenience to analysis the second ordinary differential equation 
(3.11) 

where < C < 1 a-nd t 7^ 0. In fact the real t will be taken in (0,oo). The existence of solutions to (3.11) will 
be proved by means of shooting method. Hence, the boundary condition at infinity is replaced by the condition 
^'(0) = d, where d is a real. For any d the new initial-value problem has a unique local solution 9d defined in 
its maximal interval of existence (0, Td), Td < 00. We shall see that for an appropriate d the solution 6d is global 
and satisfies 

(3.12) Odioo) = 1. 

A simple analysis in the phase plane reveals that problem (3.11) may have solutions only for r > 0. In fact the 
ordinary differential equation in (3.11) is considered as a nonlinear autonomous system in R^, with the unknown 
(6*, 9'), mainly 




(3.13) 



ip' = —Tip + 1 — 



g2 



subject to the boundary condition 

(3.14) 0(0) = C, ^(0) = d. 
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The linear part of the above system at the equihbrium point (1,0) is defined by the matrix 



J 





-2 



The eigenvalues of J are 

if T > and for < r < ^/8, 



Al — n ' ^2 — 



— T — i\/8 — — T + 'j\/8 — 
Ai = 1 ^2 = 



Therefore, the hyperbolic equilibrium point (1,0) is asymptotically stable if r is positive and unstable for r is 
negative. In particular problem (3.11) has no nontrivial solutions if t < 0. If t > we deduce from the above 
that there exists 5 > Q such that for any d and C such that + (C — 1)^ < 5"^ the local solution 9d is global 
and satisfies (3.12). In the following we construct solutions to (3.11) where the condition + (C — 1)^ < 6^ 
is not necessarily required. For a mathematical consideration the parameter C, will be taken in (—1, \/3]- The 
following theorem deals with nonnegative values of C,. 

Theorem 3.1. Let < C < and d be a real such that 
(3.15) d2<2C(l-^). 

Tiien the local solution 9^ is global, nonnegative and tends to 1 as t approaches infinity. 

The key point of the proof of this theorem is to find d such that 6d{t) > for all i > 0. To show this we consider 
again the function E{9{t), ip{t)) = |<(f (0^ + j^'"^ ^ ^- Along an orbit we have 

j^E{9{t),^{t)) = -Mtr<0. 

Hence 

Ei9dit),9',{t))<E{Cd), 

as long as 9d{t) exists; that is t < Td- The following result shows that 9d>0 on {0,Td) and then Td = oo. 

Lemma 3.1. Let < C < V3 and d satisfying (3.15). Then 9d is nonnegative, global, bounded and its Erst 
derivative also is bounded. 

Proof. First we note that from (3.13) and (3.14), there exists to > 0, small, such that 9d is positive on 
{0,to). Assume that 9d vanishes at some ti > to and suppose that 9'^{ti) ^ 0. Because 

E{C,d) > Ei9dit),9'dit)) >h'^{t^f, 

for all < t < ti. we deduce ^d^ > ({1 — ^C"^), which contradicts (3.15). Therefore 9'^{ti) = 0. In this case 
we deduce from the equation of 9d that 9'^{ti) ~ 1 and then 9d is nonnegative on a some neighborhood of ti. 
Consequently the local solution is nonnegative as long as there exists. To show that 9d is global we note that 

EiC,d)>^9'^itr + ^9l{t)~9dit)>~^, 

for all t < Td, since 9d is nonnegative. Hence 9d and (then) 9'^ are bounded. Consequently 9d is global. The 
lemma is proved I 
Now, we are ready to prove Theorem 3.1. More precisely we have. 
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Proposition 3.1. Let t > 0,C e [0,^3] and d such that (3.15) holds. Then, the global solution 9d to 
(3. 13), (3. 14) satishes the boundary condition (3.12) at inhnity. 

On account of Lemma 3.1 we need only to show that 9d goes to 1 as i approaches infinity. Since this result is 
the broad goal of the present section, we give two proofs. 
The first proof. 

Using the fact that 6d and 9'^ are bounded we deduce from the equation of 9d that 6'^ is bounded. On the 
other hand, since ■^E{6d,9'^) = —rO'j^, the function 9'^ is square integrable. Now, we use the identity 

9'd{tf = 9'd{Qf + i f 9'd{sf9'^{s)ds, 
Jo 

to show that 0^(t) has a finite limit as t tends to infinity and this limit is zero. Next, we get, by differentiation 
(3.11)i 

(3.16) C = -t9" - 29d9'd. 

Multiplying equation (3.16) by 9'^, integrating the equation obtained over (0, i), we get 

r 9'^{sfds ^ l9'^{0f + 9d{0)9'M' ^9'^it)' Odmitf + 9',{sfds. 

This implies 9'^ E L^(0,oo). Since 9'^' is bounded, by using (3.16), one sees limt^oo 6'^'(t) = 0. Finally, we 
deduce from (3.11)i that 9d{t) goes to 1 as i tends to infinity since 9d is nonnegative. 
The second proof. 

This proof uses the Bendixson Criterion. Let T be the trajectory of {9d, 9'^) in the phase plane (0, oo) x M for 
< > and let r^(T) be its w-limit set at +oo. From the boudedness of T it follows that V^{T) is a nonempty 
connected and compact subset of (0,c») x R ( see, for example p[] p 226]). Moreover (—1,0) ^ r+(T), since 9d 
is nonnegative. Note that if r"'"(T) contains the equilibrium point (1,0) then r+(T) = {(1,0)}, since (1,0) is 
asymptotically stable. Assume that (1,0) ^ r+(T), Applying the Poincare-Bendixson Theorem |14l p 44] we 
deduce that T~^{T) is a cycle, surrounding (1,0). To finish the second proof we shall prove the nonexistence 
of such a cycle. To this end we define P{9, ip) ~ ip, Q{9, ip) = —rip + 1 — 9'^ , ip — 9'^ and 9 = 9d. The function 
{9,(p) satisfies the system 9' = P{9,(p), ip' = Q{9,(p). Let D be the bounded domain of the (0, (p)-plane with 
boundary F"*". As P and Q are regular we deduce, via the Green-Riemann theorem, 

(3.17) / / {d^Q + deP) dipd9 = / {Qd9 - Pdp) = 0, 

J Jd Jr+ 

thanks to the system satisfied by {9, ip). But d^^Q + dgP — t which is positive. We get a contradiction. ■ 
In the same way as in the proof of Proposition 3.1, we can see that any global solution to (3.13) which 
is bounded from bellow by some a > — 1 tends to 1 at infinity. Therefore, to complete our analysis, we shall 
determine a domain of attraction of the critical point (1, 0). Let 

^={(C,rf) e]R2:C>-i,id2 + cQc'-i) <^}. 

Proposition 3.2. For any {d, C) in V the local solution to (3.13), (3.14) is global and converges to 1 at inhnity. 

Proof. Let us consider a one-parameter of family of curves defined by 

Ei9,^)^^^' + ^9^-9^C, 

where C is a real parameter. Note that, in the phase plane, this family is solution curves of system (3.4). The 
curve ipp' ~ 29 — ^9^ + |, corresponding to C = |, goes through the point (2,0) and has the saddle (—1,0) 
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(t = 0) as its a and ui-limit sets, see Fig 3.2. We denote this solution curve by 7i, which is, in fact, an homocUnic 
orbit and define a separatrix cycle for (3.4). We shall see that the bounded open domain with the boundary Ti 
is an attractor set for (1,0) of system (3.13) where r > 0. This domain is given by E{6,ip) = C,9 > —1, for 
all — I < C < |, which is V. As -^E < any solution, with initial data in V cannot leave V, (see the proof 
of Theorem 3.1 ). By LaSalle invariance principle we deduce that for any {C,d) in V the w-limit set, T~^{(^,d) 
is a nonempty, connected subset oi V O {(p = 0} , (see ^ p 234].) However, if6'7^1,(/3 = 0isa transversal of 
the phase-flow, so the ii)-limit set is {(1, 0)} . This means that 7^ is a domain of attraction of the critical point 
(1,0). ■ 
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Fig, 3.2 A basin of attraction V of the critical point (1, 0). 



4. Numerical results 

In this section numerical solutions of the boundary-value problem (3.11) are obtained by using the fourth- 
order Runge-Kutta integration scheme with the shooting method, velocity profiles of the dimensionless velocity 
9 are plotted in term of the similarity variable t, for various value of the shooting parameter d (the dimension- 
less skin- friction), FigA.l (1) and (2) show that the numerical results are in good agreement with the above 
theoretical predictions. 
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Fig. 1.4 Velocity profiles in terms of d = e'(0) for fixed C = 0.2 and t = 0.5 (1), and t = 2.83 (2). 



5. Conclusion 

In this work the laminar two-dimensional steady incompressible, boundary layer flow past a stretching 
surface is considered. It has been shown that the problem has solutions having a similarity form if the velocity 
distribution outside the boundary layer is proportional to x™, for some real number m. In the second part of 
this paper, we are interested in question of existence of solutions t in the case where the external velocity is an 
inverse-linear function; m = — 1. This situation occurs in the case of sink flow. To obtain exact solutions the 
stream function -0 is written under the form 

(5.1) 0(x, y) ^ ^/JlU^fit) + \og{x). 

It is shown that the ordinary differential equation satisfied by / has multiple solutions for any positive 
(suction) and no solution can exist if V^, < (injection). A sufficient condition for the existence is derived: 

(5.2) C>-i, ^■^"(0)' + c(y-i) <| 

We have obtained two family of solutions according to r = {pUoo)~^^'^ . If t > \/8, /' is monotonic and goes 
to unity at infinity, but if < r < VS, we have a stable spiral. The function /' oscillates an infinite number of 
times and goes to 1. So if we are interested in solutions to (3.9) such that 

-!</'<! 

we must take Uw, and Uoo > satisfying ~Uoo < Uw < Uoc and Vw > {^vVocY^"^ ■ 

Condition (5.2) indicates also that for the same positive value of the suction parameter the permeable wall 
stretching with velocity UwX~^, Uw > has multiple boundary-layer flows. Every flow is uniquely determined 
by the dimensionless skin friction /"(O) which can be any real number in the interval 

(-^| + 2C(i-C73), y'| + 2C(i-CV3)) , 

where ( = UwU^. The case Uoo = was considered by Magyari, Pop and Keller [23 . The authors showed, by 
numerical solutions, that the boundary layer flow exists only for a large suction parameter (r > 1.079131). 
The existence of exact solutions of the Falkner-Skan equation under the present condition was discussed by 
Rosenhad [281 pp. 244-246] who mentioned that these results may be obtained by rigorous arguments which, 
in fact, motivated the present work. We note, in passing, that it is possible to obtain solutions if the the skin 

friction satisfies |/"(0)| > + 2C (1 - CV3)- 
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